Quark and Gluon Condensates in Isospin Matter 
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Applying the Hellmann-Feynman theorem to a charged pion gas, the quark and gluon condensates at low 
isospin density are determined by precise pion properties. At intermediate density around f%rrix, from both the 
estimation for the dilute pion gas and the calculation with Nambu-Jona-Lasinio model, the quark condensate 
is strongly and monotonously suppressed, while the gluon condensate is enhanced and can be larger than its 
vacuum value. This unusual behavior of the gluon condensate is universal for Bose condensed matter of mesons. 
Our results can be tested by lattice calculations at finite isospin density. 
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While the vacuum structure of Quantum Chromodynam- 
ics (QCD) is rather complicated, quark and gluon conden- 
sates may be sufficient to determine the spectral properties of 
many hadrons. The lowest dimensional condensates in vac- 
uum are H (qq)o ^ -(250 MeV) 3 and (G 2 ) (360 MeV) 4 
with the definition G 2 = a s /nG^ u G% u , where a s is the QCD 
coupling, q the light quark field, and G" the gluon field ten- 
sor. On the other hand, the in-medium behavior of the con- 
densates are of great importance for us to understand how the 
hot and dense environment modifies the vacuum structure and 
the hadron properties. Especially, the quark condensate is the 
order parameter of chiral symmetry restoration, and the gluon 
condensate may be related to the deconfinement phase transi- 
tion n. 

For systems at zero temperature but finite baryon density 
Pb, the ratio R q (ps) = (<??)p B /(<7<z)o an d the difference 
Dg(ps) = (G 2 ) PB — (G 2 )o between the condensates in baryon 
matter and in vacuum can be expressed as [ 1] 
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where rn n , f n , on and S are pion mass, pion decay constant, 
nucleon a term and strangeness content of nucleon in vacuum, 
and ■ • • denotes the higher order correction. Taking only the 
linear terms of R q and D g which are valid at low density, the 
quark condensate at nuclear saturation density is 25 — 50% 
smaller than its vacuum value, but the gluon condensate is 
reduced by only 3 — 6%. The large uncertainty is from the 
(7 term which is not yet precisely determined. To obtain the 
behavior of the condensates at high density, effective models 
may be used. 

Recently, the study on QCD in medium is extended to finite 
isospin and strangeness densitiesyl |4J. The physical moti- 
vation to discuss isospin matter is to understand the mecha- 
nism of QCD phase transitions at finite density. While there 
is not yet precise lattice result at finite baryon density due 
to the fermion sign problem, it is in principle no problem to 
do lattice simulation at finite isospin density |3|. The QCD 
phase structure in isospin matter is also investigated in many 
low energy effective models, such as chiral perturbation the- 
ory J3,|lt], ladder QCD [6], random matrix method [7], strong 



coupling lattice QCD |8|], and Nambu-Jona-Lasinio (NJL) 
model H [M d 0]. Very recently, the isospin matter is 
discussed in the frame of AdS/CFT 111 311 . In this paper, we 
study the quark and gluon condensates in a model independent 
way at low isospin density and compare the prediction with a 
model calculation at high isospin density. While it is generally 
believed that the density effect will reduce the condensates, we 
will see that the gluon condensate at finite isospin density can 
be enhanced and even be larger than its vacuum value. Our 
result is qualitatively different from the case at finite baryon 
density, and the physical reason may be the repulsive meson- 
meson interactions. 

Since pions are the lightest excitations of QCD carrying 
isospin, the ground state of QCD at small isospin density can 
be considered as a dilute gas of charged pions. Without loss 
of generality, the matter is assumed to be composed of 7t + 
mesons with isospin density pi > 0. The relevant 7r — it 
interaction in such a system is characterized by the s-wave 
7T — 7T scattering length a in I = 2 channel. The value 
of a is predicted many years ago [ 14] and can be taken as 
a = m w /(16wf%) ~ 0.043m- 1 . 

At zero temperature and low isospin density, when the con- 
dition p\a 3 1 is satisfied, the pion matter which is in Bose- 
Einstein condensation state can be described by the LHY 
(Lee, Huang and Yang) theory 1U5I1 which is accepted as a 
general theory for weakly interacting Bose gas II 1611 . The en- 
ergy density of the system is written as 
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The first term on the right hand side corresponds to the rest 
energy, and the second one is the interacting energy in LHY. 
Note that the leading order correction 2irap 2 /m 7T was pre- 
dicted by Bogoliubov IU7I1 before LHY. The pressure V and 
isospin chemical potential p,\ can be easily obtained from £, 
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From these relations, one can determine the pion properties 
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with the equation of state of the system, 
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The behavior of the in-medium quark and gluon conden- 
sates at low density can be derived in a model independent 
way. In the QCD Hamiltonian density Hqcd, chiral sym- 
metry is explicitly broken by the current quark mass term 
W m ass = 2m q qq + m s ss + • • • , where s is the strange quark 
field, m q and m s are light and strange quark masses, and • • • 
denotes heavy quark (c,b,t) contribution which is irrelevant to 
our discussion. Possible isospin breaking effect is neglected 
here, since it will not change our result. According to the 
Hellmann-Feynman theorem of quantum mechanics iflal . we 
obtain 

2m q {yj\ J d?xqq\iP) =m fl ^-(^| J ^xHqcdW, (5) 

where \ip(m q )) is the ground state of the system as a function 
of light quark mass, and we have multiplied the equation by 
m q to obtain renormalization-group invariant quantities. 

Applying the Hellmann-Feynman theorem (O to the isospin 
matter and to vacuum and taking into account the uniformity 
of the system, we have 
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where the derivative d£ /dm q is taken at fixed density. 

The gluon condensate can be obtained by considering the 
trace of the energy-momentum tensor T£ = — 9/8G 2 + 
2m q qq + m s ss, where we considered only the u, d and s 
quarks and neglected the heavy quark contribution. The dif- 
ference between the expectation values of the trace of the 
energy-momentum tensor in isospin matter and in vacuum is 
(T£) PI - (Tfto = £-3V, which leads to the following result 
for the change in the gluon condensate, 
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where we have ignored the strangeness content of pions, 
m s ((ss) Pl - (ss) ) = 0. 

Taking the Weinberg result a oc m^/f 2 IU411 and neglect- 
ing the m g -dependence of the decay constant f„, we have 
da/drriir = a/m n . Since the energy density £ is only a 
function of m, at fixed pi, the derivative in (0 can be ex- 
pressed as d£/dm q — (d£/dm n )(dm 7V /dm q ). Combin- 
ing with the Gellmann-Oakes-Renner relation 2m q (qq)o = 
—m 2 f 2 and the LHY energy density (fJJ, we finally obtain 
the ratio R q {p\) = (qq) pj (qq)o an d difference D g {p{) = 
(G 2 ) PI — (G 2 )o between the condensates in isospin matter and 
in vacuum, 
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The linear pi-dependence of R q and D g can be completely 
determined via only two parameters, the pion mass and 
decay constant / w . From the quark condensate, it is natural 
to define a density scale po — f%m„, which is approximately 
equal to the nuclear saturation density p sat (po ~ l.lp sat ), and 
the linear dependence is valid for pi -C po. The higher or- 
der correction in isospin matter is quite different from the one 
in baryon matter. For the quark condensate, the leading order 
correction in isospin matter is 0(p^ 2 ), while it is 0(p^/ 3 ) [1] 
in baryon matter due to its fermionic nature. For the gluon 
condensate, since the coefficient m^/2 of the linear term in 
isospin matter is much less than the one (m^ — ctn — S) in 
baryon matter and the tt — it interaction is repulsive, the com- 
petition between the linear and leading terms may make the 
condensate to increase at high isospin density. Neglecting the 
term 0(p^ 2 ) which is indeed small even at high density and 
using the Weinberg result a = m n / (16tt / 2 ) [14], the gluon 
condensate starts to increase at p\ = po and becomes larger 
than its vacuum value at p\ > 2po- While this prediction is 
beyond the validity density region pi <C po of the LHY equa- 
tion of state, we do expect that the gluon condensate may go 
beyond its vacuum value at high isospin density. We will ex- 
amine it in the following with a chiral quark model. 

To estimate the behavior of the quark and gluon conden- 
sates at high density is beyond the above approach, since 
the composite nature of pions may become important at high 
isospin density. Especially, at extremely high p\, the isospin 
matter is expected to become a weakly coupled Fermi super- 
fluid JH. In this case, the element constitutes of the system are 
no longer pions but quarks. To have a complete understanding 
of the quark and gluon condensates at finite isospin density, 
we adopt an effective chiral model at quark level to describe 
the evolution from a weakly interacting Bose condensate to a 
Fermi superfluid. 

One of the models that enables us to see directly how the 
dynamic mechanism of chiral symmetry breaking and restora- 
tion operate is the NJL model IU9I1 applied to quarks B20I1 . Re- 
cently, this model is extended to finite isospin chemical po- 
tential isil [Tol [Til [l2tl . The Lagrangian density of the model 
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where ip = (u, d) is the two-flavor quark field. The current 
quark mass m q , the coupling constant G and a high momen- 
tum cutoff A due to the non-renormalization of the model are 
phenomenological parameters and can be determined by fit- 
ting the pion mass, pion decay constant and quark condensate 
in vacuum. The isospin density enters the model via intro- 
ducing an isospin chemical potential fij, corresponding to the 
isospin charge I — J d 3 x-ipjo T^ip. The order parameters for 
chiral symmetry breaking and isospin symmetry breaking are 
respectively the quark condensate (qq) and pion condensate 

At mean field level, the thermodynamic potential of the 
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isospin matter in the NJL model can be evaluated as [ 1 
d 3 k 
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M 2 , where M = m Q 



+ A 2 , £± = e±/xi/2and 
AG(qq) is the in-medium 



with the dispersions B-t 

quark mass, Mo is its vacuum value, and A = —4G(uij5d) 
is a BCS-like energy gap. The physical quark and pion con- 
densates correspond to the minimum of the thermodynamic 
potential, 
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which together with the isospin density 
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determine self-consistently M, A and pi as functions of pi. 

From the above gap and density equations, a nonzero 
isospin density is associated with a nonzero pion conden- 
sate 11211 . In the limit of pi — > 0, we have analytically 
pi — » m^. In Fig[TJ we show the numerically calculated 
isospin chemical potential pi and pressure V = — f2 in the 
NJL model and compare them to the LHY result (O with the 
standard Weinberg value a = m^KlQuf 2 ). Note that the 
higher order correction 0{^J pia 3 ) in LHY is too small to be 
observed when the isospin density is not high enough. We find 
a very good agreement of the two calculations at sufficiently 
low density. At p t ~ O.lpoi the composite nature of pions 
emerges and the NJL model starts to deviate from the LHY 
result. Note that for baryon matter, the NJL model agrees 
with the low-density result (fTJ only for some special model 
parameters. 

Having the correct low density limit, we then consider the 
quark and gluon condensates at finite isospin density. The 
quark condensate is obtained via either directly solving the 
NJL gap equations ( II \\ or the Hellmann-Feynman result (O 
through the NJL energy density £ = — V + fiipi, and the gluon 
condensate can be calculated via Eq. (0 together with the 
Hellmann-Feynman result ©. The numerical results from the 
LHY theory and NJL model are shown in Fig |2] for p\ up to 
2po- At very low density, the NJL model agrees well with 
the model independent result, but at intermediate density, the 
LHY theory for the dilute pion gas fails and the difference 
between the two results becomes significant. The quark con- 
densate decreases much faster than the gluon condensate. At 
p\ = po, the change in the gluon condensate is very small, 
but the quark condensate has already been reduced by about 
45%, which is almost the same as the condensate suppres- 
sion in baryon matter at nuclear saturation density. The other 
significant characteristic of the gluon condensate is its non- 
monotonous density dependence. While this property is out- 
side the validity density region of the LHY theory, it is con- 
firmed in the NJL model. The gluon condensate firstly drops 




FIG. 1: The isospin chemical potential pi and pressure V at low 
isospin density, calculated with the NJL model and LHY theory. 



down, then turns to go up at p\ ~ 0.6po in the NJL model 
and pi ~ po in the LHY theory, and finally becomes larger 
than its vacuum value at pi > 1.14po and pi > 2po in the two 
calculations. 

The strong suppression of the quark condensate and the en- 
hancement of gluon condensate at intermediate density will 
induce the evolution from a weakly interacting Bose conden- 
sate to a strongly interacting Fermi superfluid, i.e., the so- 
called BEC-BCS crossover 112 111 . From the fermion excitation 
spectrum, a Fermi surface is opened with pi/2 > M J22ll at 
pi > 1.6/30- This type of strongly interacting Fermi super- 
fluid may be a realization of the quarkyonic matter proposed 
recently l23ll . In such a matter, the chiral symmetry is approx- 
imately restored, but the quarks are still confined. Our pre- 
diction on the non-monotonous gluon condensate in isospin 
matter is consistent with the expectation that the hadron and 
quark phases are continued and there may exist no deconfine- 
ment phase transition at zero temperature yO . 

In summary, we have investigated the quark and gluon con- 
densates in isospin matter. At low isospin density, we derived 
the model-independent condensates by taking the Hellmann- 
Feynman theorem for a dilute pion gas. Unlike the baryon 
matter, the condensates in isospin matter are determined by 
precise pion properties. Both the low density formula and the 
NJL model calculation show that the repulsive interactions be- 
tween the charged pions may induce an unusual behavior of 
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the gluon condensate: it decreases slightly at low density and 
turns to increase at intermediate density. While there are no 
lattice data of the gluon condensate at finite isospin density, 
we expect our conclusion can be examined in the future lat- 
tice calculations. 



Our calculation for isospin matter can be easily extended 
to QCD at finite isospin density p\ and strangeness density 
ps which is associated with kaon condensation [4]. For 
kaon matter with p$ = 2pi, taking into account the relation 
m K f K — — (m s + m q )({ss)o + (qq)o)/2 for the strangeness 
condensate (ss)o in vacuum and the kaon mass mx and kaon 
decay constant fx, while the light and strange quark conden- 
sates behave differently, the gluon condensate in ([8]) for pion 
matter is valid for kaon matter, if we replace m w , p\ and a by 
itik , Ps an d kaon-kaon scattering length in 7 = 1 channel, 

D g {ps) = -8/9 (m K p s /2 - ^a K pl/m K H )■ Taking 

the recent lattice QCD result tjikclk = 0.352 l24ll . the turn- 
ing density where the gluon condensate starts to increase is 
ps ~ f K rrtK for kaon matter which is the same as p\ ~ f^m„ 
for pion matter. 



FIG. 2: The condensate ratios R q — (qq) Pl / {qq)o and R g — 
(G 2 ) P J (G 2 )o from low to intermediate isospin density, calculated 
in the NJL model and LHY theory with <G 2 } = (360 MeV) 4 . 
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